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NON-EXPLOSION CRITERIA FOR RELATIVISTIC DIFFUSIONS 

ISMAEL BAILLEUL AND JACQUES FRANCHI 

Abstract. General Loreiitz covariant operators, associated to so-called 8 (or S)-relativistic 
diffusions, and making sense in any Lorentzian manifold, were introduced by Franchi and 
^D ' Le Jan in |FLJ07| . |FLJ10| . Only a few examples have been studied. We provide in this 

^N , work non-explosion criteria for these diffusions, which can be used in generic cases. 
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1. Introduction 

It is well known that the metric completeness of a Riemannian manifold does not prevent 
Brownian motion from exploding within a finite time with positive probability. The situa- 
tion is now well-understood, in particular thanks to the works of Yau, |Yau78] . Grigor'yan 
|Gri86] . Takeda |Tak89) . |Tak91] . and very recently Hsu and Qin |HQ10| , to cite but a few 



names. Different lines of approach have been used. Yau and Grigor'yan treated the analytic 
counterpart of the completeness problem and investigated the well-posedness of the para- 
bolic Cauchy problem, the former using local information on the geometry under the form 
of curvature bounds, the latter using a global information under the form of an upper bound 
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2 I. BAILLEUL AND J. FRANCHI 

for the volume of large balls. Takeda used a purely probabilistic method due to Lyons and 
Zheng in |LZ88] . based on reversibility. This approach was recently improved by Hsu and 
Qin in |HQ10| . Hsu used stochastic analysis in [ Hsu02) . Theorem 3.5.1, to control the radial 
process, by estimating the Laplacian of the distance function to a fixed point in terms of 
curvature bounds. All these results are tied down to the metric framework provided by a 
complete Riemannian manifold. 

A natural analogue of Brownian motion in a Lorentzian setting was first introduced by 
Dudley in |Dud66j in the special relativistic case, and extended to the general relativistic 



framework by Franchi and Le Jan in |[FLJ07j. It belongs to a larger class of relativistic 
processes introduced in |BailO| and |FLJ 10J. defined in purely geometric terms, and collec- 
tively refered to as relativistic diffusions. Their trajectories represent the random motion in 
spacetime of a small massive particle, and make sense only in the unit tangent bundle or 
in the orthonormal frame bundle. Only a few examples have been studied in detail up to 
now: in Minkowski spacetime (the framework of special relativity) [Dud66] . |Bai08) . [BROS] , 
in Robertson- Walker spacetimes (models of universe with a big-bang) |Ang09| , Godel space- 
time (a causally paradoxical universe) |Fra09] . and Schwarzschild spacetime (a model for an 
isolated star or a black hole) |FLJ07] . 

Apart from the works [BailO] and [FLJIO], no general study of these intrinsic random 
processes was done. As a first step towards a better understanding of these processes and 
their interplay with the geometry of the ambient spacetime, we provide in this work some 
non-explosion criteria for some generic classes of Lorentz manifolds. In addition to being 
a natural question, the completeness issue is strongly related to important questions in 
general relativity. Indeed, dating back to Penrose and Hawking's incompleteness theorems, 
the appearance of singularities in Einstein's theory of gravitation has been recognized as 
unavoidable under quite natural assumptions. Although there is no agreement on what 
should be called a singularity of a spacetime, the existence of incomplete geodesies has been 
widely used as an indicator of such a singular feature. In so far as the random dynamics 
considered in this work ( §2.21) can be seen as intrinsic perturbations of the geodesic flow, 
their completeness/incompleteness is a distinguishing feature of a spacetime. We refer the 
reader to |Baill] for a first approach of stochastic incompleteness. 

The paths of the random processes we shall consider are (almost-)all C^ paths parametrized 
by their (proper time) arc length. What could possibly make them explode? In a complete 
Riemannian manifold, any such path would have to be at time s in a closed ball of radius 
s with centre its starting point, so it cannot explode. There are two problems with the 
Lorentzian setting: a Lorentzian manifold has no metric or finite distance function associated 
with its structure, and the set of unit tangent vectors at any point is non-compact. As a 
result, even in Minkowski spacetime, one can construct exploding paths with finite (proper 
time) arc length. 

To start our investigations, we shall take advantage in Section |3] of the bundle structure 
of the state space of the process, to exhibit a one-dimensional sub-process whose control 
is possible in the class of globally hyperbolic spacetimes. This structure allows indeed to 
define some Lyapounov function and leads to a non-explosion criterion by using a simple 
and well-known observation due to Khasminsky. 

With a metric missing, the completeness notion used in a crucial way in the Riemannian 
setting becomes unavailable. Busemann, Hawking and Ellis, Schmidt, Beem and Ehrlich, 
proposed different notions in replacement. Schmidt's idea is to give a Riemannian structure 



NON-EXPLOSION CRITERIA FOR RELATIVISTIC DIFFUSIONS 3 

to the orthonormal frame bundle. We consider Schmidt 6-completeness notion in Section HI 
showing how it leads to a stochastic completeness result for some of the relativistic diffusions. 

This result can be significantly improved by adapting Takeda' strategy |Tak91) . as im- 
proved by Hsu and Qin [HQl Oj, to the Lorentzian setting. This is however far from being 
straightforward, since we are working in a non-symmetric, non-elliptic setting, where the 
main ingredients of Takeda's method (use of symmetry and reflected Brownian motion on 
the boundary of large Riemannian balls) have no obvious Lorentzian counterpart. To over- 
come this difficulty, we use in Section [5] a sub-Riemannian structure well-adapted to our 
setting, and which will somehow play for us the role of the non-existing Lorentzian distance. 

Aknowledgements. We thank E. Trelat for his guidance in the realm of control theory, 
and A. Oancea and P. Pansu for their help in proving Lemma [T5l 
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2. Relativistic diffusions 

2.1. Basic geometrical setting. Recall Minkowski space is the product 
equipped with the metric 

9Miq,q) '■= t"^ -\x^\ \x'^\ , for any g = (t, x) G M^'"', 

where {t,x^, . . . , x'^) denote the coordinates of q in the canonical basis |eo, ei, . . . , ef^} of M^''^. 

Let (M, gf) be a smooth (1 + c?)-dimensional Lorentzian manifold (with d > 2), which we 
shall always suppose to be oriented and time-oriented. (We refer the reader to the books of 
Hawking-Ellis [HE73J and O'Neill |0'N83| for the basics on Lorentzian geometry.) Given any 
point m G M, it is usual to consider an orthonormal basis {eg, ..., e^^} of the tangent space 
T^M as an isometry e from (R^'"^, (^m) to (T^M, (7^); so, strictly speaking, ej = e(ej). The 
orthonormal frame bundle of M is just the collection 

OM = {$ = (m, e) I m G M, e an orthonormal basis of (T^M, gm)}- 

We shall write OU = |$ = (m, e) \m E U, e an orthonormal basis of T^M} for any subset 
U of M. For a small enough U and a chart x :U ^ ]R^+^ on it, we shall write Gj = e'jdr^k for 
each vector e^ of a frame e ; this decomposition provides local coordinates (x*, e^) on 014. 

Each fibre O^M is modelled on the non-compact orthogonal group 0(1, c?), which has 
four connected components. We shall be interested in dynamics leaving these components 
globally fixed. We choose to consider only one of them, specified by the requirement that eg 
should be future-oriented and that the orientation of e should be direct. We shall still denote 
the resulting frame bundle by OM, as there will be no risk of confusion. The Lorentz-Mobius 
group 500(1, d), i.e. the connected component of the unit in 0(1, d), acts properly on OM. 
This natural action induces the canonical vertical vector fields {Vij)o^i<j^d- The subgroup 
of elements in SOo{l,d) that fix eo can be identified with the rotation group SO{d), and 
generates the vector fields (l^ij)i^i<jX(i- To shorten notations we shall write Vj for Voj; it 
generates boosts, that is hyperbolic rotations in each fibre, and reads in the above local 
coordinates: 

Throughout this work, TM and OM will be endowed with the Levi-Civita connection, 
inherited from the Lorentzian pseudo-metric g. Last, we denote by i^o the vector field 
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generating the geodesic flow on OM. Denoting by Ff • the Christoffel coefficients, we have 
in the above local chart on 



(2.2) if, = eSa.,_e^'elri^.^. 

We shall denote by T^M the future-oriented unit tangent bundle over M, with generic 
element {m,m). In Minkowski spacetinie M^"'^, it is the product of M^'*^ by the hyperboloid 
U= {q = {t,x) e M^'''; g{q,q) = 1, t > O}. The bundle T^M is locally modelled on that 
product. (Consult |HE73) or |0'N83) for some background.) Denote by tti the projection 
(m, e) I— 7- (m, eg = m) from OM to T^M, and by ttq the canonical projection OM — t- M. 

2.2. Relativistic random dynamics. Relativistic diffusions model the random motion 
in spacetime of a small massive particle parametrized by its proper time, providing random 
timelike paths ; so, properly speaking, their mathematical counterpart are random trajecto- 
ries {ms,rns) in the future unit bundle T^M subject to the condition ^rrig = rfis . Yet it 
happens to be more convenient to define random dynamics in the orthonormal frame bundle 
OM as it bears more structure than T^Wl ; these diffusions on OM are constructed so as to 
have a projection on T^M which is itself a diffusion. Such a construction is remniscent of 
Malliavin-Eells-Elworthy's construction of Brownian motion on a Riemannian manifold as 
the projection of a diffusion on the orthonormal frame bundle. 

2.2.1. Dynamics in OM. Given any smooth non-negative function O : T-^M — )■ M+, identified 
to a S'0((i)-invariant function on OM by setting 0($) := 0(7ri($)), consider the following 
Stratonovich differential equation on OM : 

(2.3) od$, = Hoi^s) ds + \Yl ^M^s) Vji^s) ds + ^/Qi^ ^ Vji^s) odwi , 

where w is a rf- dimensional Brownian motion and where we understand a vector field as 
a ffist order differential operator. This equation has a unique maximal strong solution, 
defined up to its explosion time (. 

It is clear on this equation that the (ei, ..., e(i)-part of $s is irrelevant in defining the 
dynamics of (?7is,eo(s)) since G($) depends only on 7ri($) ; this is the reason why this 
diffusion on OM projects down in T^M onto a diffusion. Consult |FLJ07) . Theorem 1, or 
|FLJ10] . Theorem 3.2.1, or §3.2 of |BailO) . for the details. The diffusion in OM has generator 



(2.4) ge = Ho+'^}_^ V,{QV,). 

We shall generically call these relativistic dynamics Q- diffusions (the S-diffusions of 
|FLJ10] ). These diffusions are covariant, in the sense that any isometry of (M, (?) maps 
a 0-diffusion to a 0-diffusion (with the same B : the law is preserved, up to the starting 
point), and admit the Liouville measure as an invariant measure. The TTo-projections (on the 
base manifold M) of their trajectories are almost-surely C^ paths. A G-diffusion ($s)o<s<c 
solving Equation (12. 3 p is parametrized by proper time s > . The particular case = 
gives back the deterministic geodesic flow, and the case of a non-null constant O gives back 
the relativistic diffusion as defined ffist in |FLJ07] . which we shall call the basic relativistic 
diffusion. It is described in simple terms in Minkowski spacetime. Although the metric 
qm is non-definite positive, its restriction to any tangent space of the half sphere EI of unit 
tangent vectors is definite negative; this turns H into a Riemannian manifold with constant 
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negative curvature. Dudley's diffusion [iris, ^s) = [f^s, i^o{s), • • • , 6(^(5))), which is the basic 
relativistic diffusion in Minkowski spacetime, corresponds to taking rris = rriQ + J^ eo(r) dr, 
and for the velocity eo(r) a Brownian motion on H. The remainder ei(r), . . . , ed{r) of the ba- 
sis is obtained by paralell transport of ei(0), . . . , e^(0) along the Brownian path (eo(M))os:Msgr- 
The following elementary lemma proved in |Baill| . §2.2, gives an intuitive picture of the 
0-diffusions, for depending only on ?7i G M. 

Lemma 1. Let 7 : [0,T] -^ Wl be a C"^ timelike path parametrized by its proper time, and 
To e OM such that tiiIVq) = (7(0), 7(0)) G T^M. Then there exists a unique C^ path 
(\l/s)„^ ^^ in OM, and some unique C^ real-valued controls h^, . . . jh"^ defined on [0,T], such 
that io = To, vri(^,) = (7(5), 7(3)) and 

d 

So the 0-diffusion is obtained by replacing the deterministic controls of a typical C^ timelike 
path by Brownian controls with position dependent variance Q{ms). 

On a manifold with non-positive scalar curvature R, taking G($) = —g^R (for a non- 
null constant q), one gets a dynamics which can be truly random only in non-empty parts 
of spacetime; it was called R-diffusion in |FLJ10) . Denote by T the energy-momentum 
tensor of the spacetime. Taking 0($) = f)^T(eo,eo), we get what was named the energy 
diffusion in [FLJIO) . See [BailO] for more general models of diffusions. 



2.2.2. Dynamics in T^M. Denote by V^ the gradient on T^M, identified with the hyperbolic 
space tf by means of the metric g^ , and by Cq the vector field generating the geodesic flow 
on T-'^M . Note that Ttti{Hq) = Cq and Ttti{Vj) =: V*^ = e^, d^k (with Einstein summation 
convention). The projection on T^M of the OM- valued diffusion has the following SO{d)- 
invariant generator : 

For a constant B the operator Cq has the following expression in the local coordinates 
introduced in §2.1[ 

£0 + 1 A^' = m'^ +[lQm'- MM rj(m)) ^ + | [rn^rn' - g'\m)) ^' 
where A'' denotes the vertical Laplacian. We have for a generic G : 

(2.5) Ce = C, + %A^ + \ {m'm' - g'\m)) ^ ® ^ 



dM dM 



2 2V J Q^k Q^i 

The purpose of this work is to provide some conditions under which the 6-diffusions have 
almost-surely an infinite lifetime ( . In so far as we are mainly interested in the T^M-valued 
0-diffusions as models of physical phenomena, while we shall mainly be working with OM- 
valued diffusions, it is reassuring to have the following fact, which essentially means that the 
possible explosion of ($s)o<s«^ is never due to its (ei, ..., e(i)-part. 

PROPOSITION 2. The Q-diffusion on OM and its T^Wl-projection have the same lifetime. 
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Proof - Write $« = (nis; (m^,ei(s),. . .,ed{s))] e OM and 0^ := 7ri($s) = {ms,rns) e 

T^M.. Using the local coordinates (x*^, e^)o^A;/<d;i^jXd, Equation (12 .3^ defining the 0- 
diffusion reads: 

QQ 

dm'; = dM^ - rf^(m,) m^ m^ rfs + f 6(0^) m'; ds + ^ ("^s"^f - fl'''^("^s)) 0^-7(0^) ds , 
de';{s) = v/e(0,)mj dwi - rf,(m,) e^^(s) m*, ds + i 6(0,) ej(s) ds + \ ^,6(0,) mj ds , 



with the martingale term dM'^ := a/O(0s) ej(s) dwi- (See Section 3.2 of |FLJ10| for the 



(^]{.s)m'l gq(,{ms 



computation of the Ito correction.) Setting Cq = m and 77*" := 77" := li=„=o — li<i=ns£d ? 
and noticing that the matrix (t]*"' e^ ^ffc^) . .^ is the inverse of the matrix (e^) ^^^ , it 
follows from the above system that we have for all < fc < (i , 1 < j < c? : 

de%s) = ml T]^ elis) vK)rfMf - rf,(m,)e5(s)m^ ds + i e(0,)ej(s) ds + i VjQ (0,)mj ds 

= - e'jis) rf^(m,) m^ rfs + I 6^(5) 9(0,) ds + I VjQ{(f)s) m] ds 

do 

dm[ + r.p(m,)m*m^ t/s - f e(0,)m^ c?s - ^ [mPmf - ^*'^(m,)] —^{(j)s)ds 

So the matrix (e^(s)) and the frame- valued diffusion ($s)os:s<c satisfy a linear sto- 

chastic differential equation, conditionally on (0s)o^s<c • It is thus well defined up to the 
explosion time C, of the T^M-valued 0-diffusion. > 

This point being clarified, we shall work freely in the sequel with 0-diffusions on OM. 

3. A FIRST NON-EXPLOSION CRITERION 

We give in this section a simple non-explosion criterion, well-suited to investigate the 
behaviour of the B-diffusions in the largely used class of globally hyperbolic spacetimes. A 
Lyapounov function is introduced for that purpose, and leads to a non-explosion criterion of 
a different nature than the typical Riemannian criteria mentionned in the introduction. 

The idea is roughly the following: if we can find a function / = /($) which has compact 
level sets {/ ^ A}, and does not increase along the trajectories, then the dynamics cannot 
explode. This was noted first by Khasminsky in a stochastic context; we state his observation 
here for the relativistic diffusions. 

Lemma 3 (Khasminsky). // there exists a non-negative function f on OM and a positive 
constant C such that Q^f ^ C f , and f goes to infinity along any timelike path leaving any 
compact in a finite time, then the Q-diffusion has almost- surely an infinite lifetime. 

Proof - The condition ^0/ ^ C f implies that the real-valued process (e^'"'*/($s)) ^ is 
a non-negative supermartingale. Denote by r„ the (possibly infinite) exit time from the 
level set {/ ^ n}. By optional stopping, we have 

/($o) ^E[e-^^"/(<l>,J] =nE[e-^^"]. 

This implies that Tn goes to infinity as n goes to infinity; as C, = lim r„ , this proves 

Khasminsky 's statement. t> 

As 0-diffusions have no a priori reason not to explode, such a Lyapounov function will 
generally not exist. Yet, it is possible to construct such a function in some classes of space- 
times of interest for cosmology and theoretical physics. We give below two such examples. 
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The construction of the function / uses the same recipe in both cases : if there exists an 
intrinsic distinguished future-directed timehke C^ vector field U G T^M, we can define 

(3.1) f{^):=gm{U^,m)- 

recall that 7ri($) = (171,771) G T^M. For this choice of /($), which is the hyperbolic angle 
between U and rn , we have / ^ 1 , and 

(3.2) Hof{<^) = Vrn{9{U,7h)) = g{V^U,77l) . 

The following lemma shows why / is a good choice to apply Khasminsky's criterion. 

d 

Lemma 4. We have on OM: | E ^i(0 ^j/) = f © / + K/"^'' - f^'') |^ ■ 

j=i 

Proof - Choose local coordinates for which U = d^o, so /($) = m° = Cg . Using (12. ip . we 
have thus locally : 

^^^-(''ss^^«4>°^^'°' ^'f-^'-f' 

and t(y,e)(V,f) =±4e^ ^ = (A"™' - /') |^ = (/ *' " C) ^ ■ 



It follows from (ED and (E^D that 



Qef = 9{VrnU,7h) + ^Qf + Ufm' -U') ^® 



2 ^ 8771'' 

Khasminsky's criterion will thus guarantee the non-explosion of the 0-diffusion provided / 
explodes along exploding trajectories and there exists a positive constant C such that 

(99 

(3.3) g{VrnU,7h) + K/m'^ -U')^-j^^{C-lQ) g{U,7h). 

In order to turn this criterion into an effective tool, we first restrict ourselves to the following 
general class of spacetimes. This inequality become s particularly simple when 6 depends 
only on the base point m G M. 

3.1. Globally hyperbolic spacetimes. This class of cosmological models is character- 
ized by the existence of a global time function (that is a function r : M — )■ R, with timelike 
gradient) such that it has connected spacelike level sets {r = t} of r, and each integral curve 
of the vector field Vr meets each leevl set of r in exactly one point. Thus M is diffeomor- 
phic to the product / x S* of an interval / and a d-dimensional manifold 5*. Without loss 
of generality, we can suppose the interval I unbonded from above. With the example of 
Minkowski spacetime in mind, we see that a given spacetime may have an infinity of time 
functions; they are not intrinsically associated with the geometry. 

Yet, we can take for vector field U in this setting the gradient of the time function r : 
m = (t, a;) G / X 5* t-T- t, so 

/($) = g{U,7h) = V^r = m° = t > 0. 

There is no hope, though, to prove Inequality (13. 3p without specifying further the model, as 
the time function is not intrinsically defined. To proceed further, we shall look at the sub- 
class of generalized warped product spacetimes , in which the time function is supplied 
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by the model and can be seen as an absolute time. These universes are globally hyperbolic 
spacetimes M. = I x S whose metric tensor has the form 

(3.4) 5fm(m,m) = a^ |m°| - hm{m^ ,rn^) , 

where rnP is the image of m G T^M by the differential of the first projection I x S ^ I 
and TTi'^ the image of m by the differential of the second projection I x S ^ S. Write 
m = (t, x) G / X 5. The function a is a positive C^ function on M , assumed to be 
bounded on any subset I' x S where /' is bounded from above, and hm is a positive- 
definite scalar product on T^S, depending on m in a C^ way. This class of spacetimes 
contains all Robertson- Walker spacetimes (hence in particular de Sitter and Einstein-de 
Sitter spacetimes, and the universal covering of the anti-de Sitter spacetime). 

Theorem 5. Let {M,g) be a generalized warped product spacetime. If the function 
T'M ^ (,n.,n) ^ V„loga - |e(™.,h) - l(,h'|| - -^^ ||) 

is bounded below, then the Q-diffusion almost- surely does not explode. 

Proof - • We first check that if the 0-diffusion has a finite lifetime C^ then /($s) ex- 
plodes at time C,~ . To that end, consider a timelike trajectory 7 = (ms,rhs)o<s<T = 
((ts,Xs),ms) in T^M, defined on some semi-open interval [0,r), and such that 

■^rris = ths and /(7s) = is is bounded above by some positive constant C. It follows 
that to < ts <tQ + CT, and hmX^siis) ^ C'^a^^ is bounded above by a constant since 
a is bounded above on (inf /, to + CT] x S. This entails that {xs)o<s<t cannot exit a 
bounded region of S, and so that 7 must be trapped in a finite union of sets of the form 
J+(r7io) l~l J~{Qj) 1 foi' some Qj G M. A such a union of sets is compact in a hyperbolic 
spacetime (see for instance |HE73) . Section 6.6), 7 is trapped in a compact set. Would 
7 explode, it would have a cluster point at which the strong causality would fail, leading 
to a contradiction as globally hyperbolic spacetimes are strongly causal ( |HE73] . Section 
6.6). 

• The condition of the theorem is a rephrasing of the local condition (13. 3p . To see that, let 
us work in a neighbourhood V = [^1,^2] x \/ of a given point m^ , and choose coordinates 
x^ on V ] this provides coordinates (t, x*) on V, which induce coordinates on T^V: for 
m eV and m G T^M, write m = rhPdt + X] rn^d^j. 

l<j<d 

Note first that since U = a^'^dt , we have 

Using Christoffel' symbols r*^^. we have 

[V^dtY = V^(a-2) 5^ + a-^m-^ r^o, 
for a G {0, ..., d} and a summation over c in {0, ..., c?}; so 

Hof = giVrnU, m) = V^^(loga-2) m" + a-^m' T^q g^^ m^ . 
The explicit formulas for the Christoffel symbols, in terms of the metric, are 

r0o = 9i(loga), T% = dA^oga), r,, = '^h''d.Aa'), T\, = \W'dth,k , 
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for i,k G {1, ...,d} and a somniation over 1 < i < d . We thus have, after simphfications, 

Hof = -2V^(foga)m° + |mY5t(log«) - ^rfi'' dt{hk)m' 

= — \m^\'^ dt log a — 2 m° m^ d^k log a — m'^ dt{hek) rnf' ■ 

Using the unit pseudo-norm relation a^ |m°p — hik{m) m^ rr/ = 1, the above equality 
becomes : 

i/o/ = -|mY dt log a - 2 mOm'^ 9^^ loga \m°\^ dtia^) , 

that is, HqJ' = —2 m° Vm log a . The statement of the theorem follows from (13. 3p . > 
This result takes a particularly simple form in the case where B depends only on the base 
point m , as is the case of the /^-diffusion. 

Corollary 6. Let M = I x S denote a generalized warped product spacetime and 
be a bounded non-negative function on M. Then the Q-diffusion does not explode if Va is 
everywhere non-spacelike and future- directed. 

Proof - The condition of Theorem |5] reads in that case: "T^M 3 {m,fn) i— )■ Vmloga 
is bounded below'''. To rephrase this condition into the more synthetic condition of the 
statement, let us work in local coordinates, (t,x) and (t, x) for m and m respectively. 

We have t = a~^chr and x = (shr)cr , for some r G M and a G T^S with |cr|/i(m) = 1. 

Define u := dfloga and v := d^loga G T^S = W^. Then the condition of Theorem [S] 
reads : wa^-'^chr — (fjcr*)shr ^ C, for any r and a. Letting r — )■ ±oo, gives a~^u > 
|fjcr*| > 0. As the constant C can be taken negative without loss of generality, the 
reciprocal is clear. Now, since max |ficr*| = |f|/t-i(m) , the condition reads: 

klh(m)=l 

a~^u > |f |h-i(m) • Finally, as V = (a~'^dt , —h^^d^i^ the vector Vloga = {a~'^u , —k^^Vj) 
has pseudo-norm (7 (V log a, V log a) = a^'^v? — IvW-ii^-. > 0. > 

This criterion applies in particular to G-diffusions in Robert son- Walker spacetimes, recover- 
ing the results of Angst, [Ang09j , who proceeded by direct analysis of the stochastic differ- 
ential equations of the dynamics. 

3.2. Perfect fluids. Our second class of examples where to apply Lyapounov's method 
to prove non-explosion will be the set of spacetimes with normal matter whose energy- 
momentum tensor T is that of a perfect fluid. They are characterized by the datum of 
a timelike vector field [/, the four velocity of the fluid, and two functions p and p on M, 
respectively the energy density and pressure of the fluid. See [HE73J . [BEE96J . We have 
then T = pU ® U + p{^g + U ® U') , or in local coordinates, 

Tij = {p + p) UiUj + p gij . 

Such a spacetime is said to be of perfect fluid type. Notice that contrarily to the globally 
hyperbolic spacetimes no topological assumption is made on a perfect fluid type spacetime. 

Godel's universe is such a spacetime. This is the manifold R^ with the metric ds"^ = 
dt^ — dx^ + \e^^^^ dy'^ — dz^ — 2e^^^dtdy, where a; > is a constant. It is a solution 
to Einstein's equation with cosmological constant cu^ and represents a pressure free perfect 
fluid. It has energy-momentum tensor T = f/ (g) f/ , where ([/,) = (-\/2tJ, 0, ■\/2i^e^^'^^, 0) 
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represents the four-velocity covector of the matter, and u is the vorticity of this field. This 
spacetime has constant scalar curvature 2u'^. See Section 2.4 in |Fra09) . As above, the 
function / is defined by Formula ( 13. ip and can be used as a Lyapounov function under 
some conditions. The computations made in Section 13.11 work equally well in that setting 
and lead to the following results. 

Proposition 7. Let (M, g) be a Lorentzian manifold of perfect fluid type, and f defined 
by Formula ^3.1^ . Suppose f goes almost-surely to infinity along any exploding timelike path. 
If there exists a constant C such that 

Hof + iQf + \{fm^ -U')^—^^C f , 
then the Q- diffusion has almost-surely an infinite lifetime. 

In the particular case of Godel universe, the gradient VU of the velocity vanishes (since 
[/' = (5o), so that H^f = 0, by Formula (13. 2^ : and / is the square root of the energy. 

Corollary 8. Let us work in Godel universe and suppose that 3G + [ff^^^^ — t^t) ^■^ 
bounded above in T^Wi. Then the Q-diffusion has almost-surely an infinite lifetime. 
This condition holds in particular if 0($) = Q{m) depends only on the base point and is 
hounded, as this is the case for the basic relativistic diffusion and the R-diffusion in Godel 
universe. 

Note that this criterion does not apply to the energy diffusion in Godel's universe. Indeed 
one can see in that case (see Section 2.4 of |Fra09] ) that the above quantity is equal to 
50 — Auj"^ and that the energy is unbounded along the trajectories of the energy diffusion. 

Remark 9. In Einstein-de Sitter spacetime the energy diffusion explodes with positive 
probability, as proved in Proposition 5.4.2 of |FLJ10| . (This Robertson- Walker universe is 
both a warped product and a perfect fluid type spacetime.) Consult [Baillj for a first study 
of stochastic incompleteness for relativistic diffusions. 

4. 6-COMPLETENESS 

The study of dynamics in the orthonormal frame bundle is not new in general relativ- 
ity and essentially dates back to Cartan's moving frame method. However, B.C. Schmidt 
|Sch71] was the first to notice that the geometry of OM itself may be used to provide a 
conceptual framework in which studying the nature of spacetime singularities. For that 
purpose, he introduced on the parallelizable manifold OM a Riemannian metric, turning 
{i^o, •••, -f^d 5 (^j)o^i<i^(i} into a Riemannian orthonormal basis, and called it the bundle 
metric, or h-metric. The completeness of this metric structure on OM can essentially be 
phrased in terms of M- valued paths. To state that fact, recall that one can associate to any 
M-valued C^ path 7 : [0,T[-> M and e G O^^^ a unique horizontal lift 7"^ : [0,T) -^ OM 
of 7, starting from (70, e), and charactarized by the properties 

— 7J e span(ifo, •••, H^), and 710(7!) = 7s, for all s E [0, T) . 

The Se-length of 7 is defined as the Riemannian length of its horizontal lift 7^; it depends 
on e e 0-y(,M. In other words, given e G 0^„M, seen as orthonormal in the Euclidean sense, 
the 5'e-length of the M-valued C^ path 7 is the Euclidean length of its anti-development in 
(Ty„M, e). Although this length depends on e, its finiteness is independent of it; we can 
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thus talk of finite S-length of a C^ path without mentioning the frame e. Note that in a 
Riemannian setting the Sg-lenth of a C^ path is its usual Riemannian length. 

Theorem 10 (Schmidt, |Sch71] ). OM is complete for the above b-metric if and only if any 
C^ path 7: [0,T) — 7- M with a bounded S-length converges in M at time T~ . 

The above completeness hypothesis is usually called b-completeness . The Riemannian 
version of this statement is trivial as the orthonormal frame bundle with its b-metric is 
complete iff the Riemannian manifold is complete. The Lorentzian situation is more involved 
as there exists (timelike, spacelike and lightlike) complete Lorentzian manifolds M which have 
an incomplete path of bounded acceleration, so OM is not b-complete, see e.g. |Ger68] and 
|Bee76] . The non-compactness of 6*00(1, d) lies at the core of this phenomenon. 

However, the Riemannian view of a Lorentzian manifold provided by Schmidt's metric 
offers a bridge to investigate some features of the latter using the tools of Riemannian 
geometry, as the following proposition shows. 

PROPOSITION 11. Let Q be a bounded function on M. Then the Q-diffusion does not explode 
if OM. is b-complete. 

One should not be confused about that statement. It does not mean that the Riemannian 
completeness of OM implies the completeness of its Brownian trajectories, which is false. 
One cannot assign an Sg-length to a Brownian path in OM as it is not regular enough. 

Proof - • Given a horizontal C-'^-path {ps)o^s<T in OM, write 7 for its projection iToop in M, 
so p = 7^. For < s < T, denote by t^^^ the parallel transport operator along the curve 
{jr)o^r^s, with inverse tJ^^ . Also, denote by {ps)o<:s<T the anti-development of 7: this 

T^gM- valued C^-path is defined for all s G [0,T[ by the formula: Ps = TQ^^.'jj.dr. 

Jo 
Last, we shall denote by pi the coordinates of pr in the frame po, and by || . ||p^ the 

Euclidean norm in (T^^M, p^)- We have by construction dps = X] Hj{ps)pids and 

0<j<d 

7s = "Tq^sPs, as well as the identity ||7s||p^ = ||Ps||po ~ X] {Ps) ■ The b-completeness 

0<j<d 

assumption means that 7 has a limit '-/j- in M at time T if 

(4.1) / \\p s\\ po ds < 00. 

Jo 

• The basic relativistic diffusion (771^,6^) is by construction the development in M 
of the relativistic Dudley diffusion in Minkowski spacetime, identified with T^^M , see 
Theorem 3.2 in [F-LJ-1]. As trajectories of the latter over a time a bounded time interval 
have almost-surely a finite length in the Eulidean norm associated with any frame of M^'*^, 
the b-completeness of OM ensures the non-explosion of the basic relativistic diffusion. 

• For a generic 0-diffusion, Formula (12. 5p implies the existence for each s G [0,C[ of an 
orthonormal basis ((y9i(s), ...,y9d(s)) of pj- in IS}'^ such that one has 

d 

dp^s = Y^ Ve(m,) ^]{s) dwi + f e(m,) p'; ds 
i=i 
for some (i-dimensional Brownian motion w. We have used the fact that B depends only 
on m to simplify the general expression. The path {ps-,Ps)o^s<c appears then as a time 
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change of Dudley's diffusion, by means of the map s i— )■ infjw I J^ Q(mj.) dr > s}. The 
result follows for a bounded function B. > 

This result can be improved in two ways: by relaxing the boundedness hypothesis on 
and by relaxing the geometric completeness assumption. The next section explains how this 
can be done in a sub-Riemannian framework by using ideas from the theory of reversible 
Markov processes. 

5. A VOLUME GROWTH NON-EXPLOSION CRITERION 

We prove in this section a non-explosion criterion involving only the volume growth of 
some sub-Riemannian boxes in OM and the function 0, as described in theorem [T^ below. 
This result is proved in Section [23] following Takeda's method, as improved recently by Hsu 
and Qin in |HQ10| . Yet, there is a real difficulty in doing this, as we are working with a non- 
symmetric, hypoelliptic diffusion, and on a principal bundle with non-compact fibres. To 
overcome these difficulties, we introduce a sub-Riemannian structure on OM, well-adapted to 
our setting, and which will somehow play for us the role of the missing Lorentzian distance. 

5.1. Sub-Riemannian framework and main results. 

5.1.1. Sub-Riemannian distance function. We have seen in §l]that the completeness of the 
natural Riemannian metric of the parallelizable manifold OM implies the stochastic com- 
pleteness of all the B-diffusions with a bounded G. One can significantly improve that 
conclusion by working with the sub-Riemannian structure on OM induced by the field of 
{d + l)-planes generated by the vector fields Hq, Vi, . . . , V^. In that setting, one can assign a 
length only to C^ paths p : [0, T] — > OM whose tangent vector belong at any time s to the vec- 
tor space spanned by Hq, Vi,...,Vd in T^^OM, say p, = p°^Ho{ps) + pl Vi(Ps)H Kpf Vd(p,). 

Such a path is said to be admissible; its length is then defined as J (X]j=o(/^s) ) ^ '^^- '^^^ 
sub-Riemannian distance between two points of OM is defined as the infimum of the length 
of the admissible paths joining these two points, with the convention inf = -|-oo. Chow's 
theorem |Cho39| ensures that the sub-Riemannian distance function T>{-,-) is finite and con- 
tinuous in its two arguments if (see e.g. |Mon0 2]) the Lie algebra generated by Hq, Vi, . . . ,Vd 
has full dimension , which holds here. Fix a reference point $i.ef ^ OM. 

(H) Completeness hypothesis. The closed boxes B\ := {V{^ref,-) ^ M o,f^ compact for 
any A > 0. 



This completeness hypothesis rules out the pathological examples of Geroch j Ger68j and 
Beem |Bee76j : it does not depend on the arbitrary choice of $o- Unlike its Riemannian ana- 
logue, the sub-Riemannian distance function P($ref, ■) is not smooth in any neighbourhood 
of $ref, |Mon02] : however, it is a viscosity solution of the equation 

\HqV\^ + \V{D\^ + ■■■ + \VdV\^ = 1 

on OM\{$i.ef} (see e.g. theorem 2 in |Dra07| : we do not use that fact in the sequel). We 
shall use that quantitative information in §5.41 under the classical form given in the following 
proposition. 

Proposition 12. Fix X > 0. One can associate to any positive constant rj a smooth function 
F : OM -^ ]R+ such that 

max|F(<l>) -I?(<l>re/, $)| ^Tj 
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and we have on Bx 

\HoF\^ + \ViF\^ + ■ ■ ■ + \VdF\^ ^ 2. 

Proof - Let us introduce the Riemannian metric g^ on OM for which Hq, Hi, . . . , H^ and 
the iyij)Q^i<j^d are orthogonal, with Hq and the Vqj{= Vj) of norm 1 and the other vectors 
of norm e~^ . Denote by T'e(-) = Pe($ref; ■) the distance function associated with g^. It 
is a 1-Lipschitz-continuous function (with respect to the distance function V^) which is 
differentiable almost-e very where, by Rademacher's theorem, and has a gradient of norm 
1 almost-every where: 

(5.1) \H^V,\'' + \V{D,\'' + --- + \VdV,\^ + e'^(Y.\H,V,\''+ J^ \V,pA=l. 

(Indeed, the set of conjugate points to $o in B\ is closed and has null measure. In 
the complementary open set the distance is attained along a unique geodesic whose unit 
tangent vector at the final point is the gradient of the distance function to $o-) The 
function D^ is easily seen to converge uniformly to ©(^ref, ") on the compact box B\ (this 
is where we need these boxes to be compact); see e.g. §§0.8. A and 1.4.D of Gromov's 
article [Gro96j . As we have almost-everywhere 

by (15. ip . a standard regularization procedure yields the conclusion. I> 

5.1.2. Main results. We use the natural volume measure on OM associated with the Lorentzian 
structure. It is defined by the formula 

VoL(t/$) = YoLM{dm) ® VoLm((ie), $ = (m, e), 

where Volm(c??«) is the Lorentzian volume measure and VoLm{de) is the image of a given 
Haar measure on 5*00(1, c?) by the identification of the fibre TCQ^(m) with SOo{l,d) (see 
e.g. [HE73 J. Section 2.8, for the Lorentzian volume measure). The volume measure Vol on 
OM is uniquely defined up to a multiplicative constant. In order to avoid some unpleasant 
pathologies, we shall make the following rather mild assumption on the causal structure of 
spacetime. 

Hypothesis. (M, g) is strongly causal. 

It means that any point of M has arbitrarily small neighbourhoods which no non-spacelike 
path intersects more than once; see |HE73] . p. 192, or |BEE96| . 

Theorem 13. Let (M, g) be a strongly causal Lorentzian manifold satisfying the complete- 
ness hypothesis (H). Set 0^ := sup 0($), for any r > 0, and suppose 

, , /"^ rdr 

5.2 / -. rT=00. 

^ ' J e, log(e,VoL(5,)) 

Then the Q-diffusion has almost-surely an infinite lifetime, from any starting point. 

Condition (15. 2p has the form of the classical non-explosion condition for Brownian motion: 
-, — r = oo, first proved by Grigor'yan in |Gri86| and has precisely that form for 

logVOL(5,,) 

bounded. Note that no topological assumption on M is needed, contrary to the results of 
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§3.11 One can give a quantitative version of the above theorem by providing an upper rate 
function. 

Corollary 14. Let Wl be a strongly causal Lorentzian manifold satisfying the completeness 
hypothesis (H). Set h{p) =pifQ = 0; otherwise, pick a constant Rq such that O/j,, > and 
set for p > 

{ f^ rdr } 

^ ' I Jr, e, iog[e,voL(i?,)] ^j 

Then, given any $o ^ OM, there exist Rq > and a positive constant C such that we have 
F,^g- almost- surely 

I)($o,$.)^C/i(Cs). 

We prove Theorem [T3] following Takeda's method, explained in the next section. To adapt 
it to our setting, we shall introduce in §5.31 a modified B-diffusion on some compact space; 
it is used crucially in the proof of Theorem [13] given in §5.41 

5.2. Takeda's method. 

5.2.1. The main ingredients. Using an idea of Lyons and Zheng, |LZ88| , Takeda devised in 
|Tak89j . [TakQl J. a remarkably simple and sharp non-explosion criterion for Brownian motion 
on a Riemannian manifold V. Loosely speaking, his reasonning works as follows. Suppose 
we have a diffusion {Xs)s^o on V which is symmetric (with respect to the Riemannian volume 
measure Vol, say) and conservative; denote by L its generator, and let / be a sufficiently 
smooth function. Denote by Pvol the measure / P^; VoL{dx) on the path space, where P^^ is 
the law of the diffusion started from x. Fix a time T > 0. As the reversed process {xt-s)o<:s^t 
is an L-diffusion under Pvol ; cipplying Ito's formula to both f{xs) and f{xT-s) provides 
two martingales M and M (with respect to the two different filtrations cr^Xg ; ^ s ^ T) 
and a{xT-s ; ^ s ^ T) respectively) such that : 



s 



f{xs) = f{xo) + Ms+ Lf{xr)dr, 

Jo 

pT-s 

f{xs) = f{xT-(T-s)) = fixr) + Mt-s + / Lf{xT-r) dr. 

Jo 

It follows that f{xs) = ^^"^"^y^^'^^ + Mi+Ml^ + J^ Lf{xs) dr, and consequently, 

f{xT)-f{xo) = '^{MT-MT). 

If ^^J" and ^^J" are bounded above, by 1 say, the previous identity provides a control of 
(/(xt) — fixo)) by the supremum of the absolute value of a Brownian motion over the time 
interval [0,T]. 

Back to the non-explosion problem for Brownian motion on V, fix a point m G V and a 
radius R > 1, and consider the Brownian motion {xs)s^o refiected on the boundary of the 
Riemannian ball B[m] R), started under its invariant measure lB{m;_R)VoL. It is a symmetric 
conservative diffusion ; denote by fB(m;R) its law. Using the Dirichlet forms approach to sym- 
metric diffusions one can apply the above reasonning to the (non-smooth, but 1-Lipschitz) 
Riemannian distance function d{m, .), which gives the estimate 



^B{m;R)(xo G B{m] 1), supd{m,Xs) = Rj ^ Yoh[B{m;R)) x 2P(sup \Bs\ > R 
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But as the Brownian motion on V behaves in the ball B{m] R) as the Brownian motion 
reflected on the boundary of B{m;R), the above inequality also gives an upper bound for 
the probability that the Brownian motion on V, started uniformly from S(?7i; 1), exits the 
ball B{m; R) before time T. Combining this estimate with the Borel-Cantelli lemma, Takeda 
proved that the Brownian motion on V is conservative provided 

liminf i?"^ log Vol (5 (m;i?)) < cx) , 

re-proving in a simple way a criterion due to Karp and Li. Takeda's method has been refined 
by several authors, culminating with Hsu and Qin's recent work [HQlOj , who give an elegant 
and simple proof of a sharp non-explosion criterion, due to Grigor'yan |Gri86| . for Brownian 
motion on a Riemannian manifold in terms of volume growth, as well as an escape rate 
function. We shall follow their method to deal with relativistic diffusions. 

5.2.2. The difficulties. The main difficulty in implementing this approach is in finding what 
can play the role of the pair "Riemannian distance function - reflected Brownian motion" 
in our Lorentzian, hypoelliptic framework. We describe in the remainder of this section a 
non-standard refiection mechanism for a Brownian motion in a Riemannian manifold which 
will serve us as a guide in the construction of the ©-diffusion reflected on the boundary of 
the sub-Riemannian boxes, as described in section 15.31 

Brownian motion reflected on the boundary of a ball B{m;R) is the simplest diffusion 
process which coincides with Brownian motion on the ball B{m]R) and has a state space 
with flnite volume. One cannot take a smaller state space if the former property is to be 
satisfled. Yet, one can make different choices if one is ready to loose the minimality property. 
To explain that fact, let us suppose that (V, (7) is a Cartan-Hadamard manifold. Given a 
point m e V let us use the exponential map exp^ at m as a global chart on V; this identifles 
the geodesic ball B{m; R) on M to the (Euclidean-shaped) ball B'{0; R) in T^V. Given e > 0, 
let us modify the metric on B'{0;R + e)'\B'{0;R) so as to interpolate smoothly between 
exp'^g on B'{0]R) and the constant metric Qm outside B'{0]R + e) (primed balls refer to 
the pull-back metric exp^ (7). Denote by g the restriction to B'{0] R + 2e) of this modifled 
metric, and deflne the compact space K as the quotient of the closed ball B (0; R + 2e) by the 
identiflcation of m' G dB (0; R + 2e) and —m'. Then the gf-Brownian motion on K coincides 
with the expj!„ (jf-Brownian motion on B'{0]R) and has a state space with flnite 5^-volume 
VoL^(K) = (1 + o(e)) YoLg(^B{m; R)). The construction of a modifled 0-diffusion given in 
section 15.31 will be reminiscent of the preceding non-standard reflected Brownian motion. 

5.3. A modified process. We start our construction of the "reflected" B-diffusion by 
constructing the compact space on which it is going to live. Fix for that purpose a reference 
point $ref e OM, the centre of the boxes Bx, and set V{^) = ©(^ref, *) for all <l> G OM. 
Fix also two positive constants A and e and consider the relatively compact open region 

U:={X<V<X + e} = Bx+e \ B^- 

Lemma 15. There exists inU a smooth hypersurface V o/OM separating dB\ from 9-Ba+e 
such that the subset Vq := {$ G V^ | Hq{^) G T$V^} is a smooth hypersurface ofV. 

The separation property means that dBx U dBx+e does not intersect V but any continuous 
path from dBx to dBx+e hits V. We thank A. Oancea and P. Pansu for their help in proving 
this statement. 
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Proof - Let us use the function F of Proposition [121 with 7] < e/A and i? > A + e , and 
fix some constants rj < ei < 62 < e/2 — r] such that Bx C {^i < F — X < 62} C Bx+e/2 ■ 
The set of regular values of {F — A) is dense in the interval {61,62) , by Sard's theorem. 
Fix a regular value c G (ei,e2), so the level set 5* := {F = c} is a smooth hypersurface 
separating dBx from dBx+e/2 ■ 

We shall now be working in W = 5* x [O, |), where we are going to construct the sep- 
arating hypersurface V as the graph of some function / : 5 — )■ [O, |), resorting to the 
transversality lemma. Denote by Gr{TU') the Grassmannian bundle over W made up 
of all the hyperplanes of TW, and associate to any function / : S* — > (O, |) the func- 
tion Gf : S ^ Gr{TU') defined by Gf{m) := {(cr, t//m(cr)) | o" G TmS}. Let "H denote 
the smooth hypersurface of Gr{TU'), made up of all hyperplanes containing Hq . Then 
Gj^{7i) is a smooth hypersurface of Graph(/) as soon as Gf is transverse to Ti . Therefore 
the statement reduces to finding a function / such that Gf be transverse to "H . 

k 

Consider for that purpose a smooth partition of unity: 15=^ dj , with {aj > 0} = 

i=i 
ipjiB") diffeomorphic under ipj to the unit ball B" C W (with u = dim(OM) - 1 = 

{d + 3)d/2j. Denoting by A the space of (the restictions to B'^ of) affine functions on W^, 
consider the map F : A"" x 5 — t- Gr{TW) defined by the formula 

G{ipi,...,ipk,m) := Gf{m), 

k 

where / = ^ a^ v^j ° ^7^- This is easily seen to be a submersion. It follows from the 

i=i 
transversality lemma that such a Gj is transversal to "H for almost-every {ipi, . . . , Lpk) G 
A^. The graph of the function / corresponding to a small multiple of such a A;-tuple has 
the properties of the statement. I> 

Let O be the set of points of the box Bx+e of the form 7(1) for some continuous path 
7 : [0, 1] — )■ Bx+e starting from a point of Bx and not hitting V] this is an open set with 
\^ as a boundary. Denote also by W another smooth hypersurface, separating V from dBx 
and transverse to Hq except on a relative hypersurface. Let now denote by O'M a disjoint 
copy of the set of past-directed frames: 

{(m, e) G GLM \ e = (eo,ei, . . . , e^) such that (m, (— eo,ei, . . . ,6^)) G OM}, 

and let O', V , Vq and W be the subsets of O'M corresponding to O, V, Vq and W. The 
equivalence relation 

(m, (eo, ei, . . . , e^)) G l^ ~ (m, (-Gq, ei, . . . , e^)) G V 

defines a manifold structure on the quotient space {OUV)VA {O' U V) j ~ , which we denote 
by £. Note that £ is compact and that its volume is in between 2 VoL(i?A) and 2 VoL(i?A+e). 
Write V for the image in £ oi V, and Vo for the image in £^ of Vq ; define the primed sets 
V and Vq accordingly. 

Remark 16. The geodesic flow is naturally well defined on £ \Vq , getting instantly from 
O to O' or from O' to O at its crossings of V \Vo . Indeed by the above definition, for 
any $ G V\Vo, either Hq{^) points outwards seen from O and inwards seen from O', or 
Ho{^) points inwards seen from O and outwards seen from O'. There is however no a priori 
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convenient way to extend the geodesic flow on Vq- This is the reason why we need to take 
care of this exceptional set. 

We define the modified relativistic diffusion on the compact manifold S as follows. 

Let a : -Ba+e — )■ [0, 1] be a smooth function equal to 1 on i?A , and whose vanishing set 
is exactly the closed part €. oi U in between W and V (this means that (t is the union 
of the trajectories (7s)se(o,i) C W of continuous paths 7 such that 70 G W, 71 G V, 
and (7s)se(o,i) does not intersect the oriented hypersurface ly U K). We extend to S the 
restiction of a to O U l^, by setting a(e') = a(e) for e' = (m, (—eg, ei, . . . , e^)) G O'M and 
e = (m, (eg, ei, . . . , ed)) G OM. We define the generator of the modified diffusion to be the 
following variant of Qq : 



d 



(5.3) g:=Ho + lJ2^,iaeVj). 

Denote by Vol^- (resp. VoLy , Vol^^/) the natural volume element on S (resp. V, W). 

Lemma 17. For Vol^- almost all starting point ^o E S, the modified relativistic diffusion is 
a well-defined S-valued process having an almost-surely infinite lifetime. 

Proof - This modified diffusion has generator Qq in Bx and in its mirror copy B'^ , and 
reduces to the geodesic flow in the region {a = 0} in between W and W . After remark 
[T6| we need first make sure that the set Vq U Vq of bad points is polar. 

Let M and M' be the orbits in the region {a = 0} of Vq and Vq by the geodesic flow. They 
have, as a consequence of lemma [151 iiuU VoL£:-measure. But as the modified diffusion 
started from any <^q & {a > 0} is hypoelliptic, its hitting distribution oiWU W has a 
density with respect to VoL^yuW- It follows that the modified diffusion, started from any 
point of ^Q\{f/ U N"'), will almost surely never hit A/" U A/"', proving that this £^-valued 
process is well-defined. 

It can behave in two ways as it approaches its lifetime: either crossing infinitely many 
times V, or remaining eventually in a compact subset of O or O'. In the latter case, its 
projection on M is a (future or past-directed) timelike path confined in a compact subset 
of O. As such it has a cluster point at which the strong causality condition cannot hold, 
preventing M from being strongly causal, a contradiction. 

In the former case, either the path eventually remains in the region {a = 0}, or it performs 
before some finite proper time an infinite number of crossings from W U W to V. Since 
the geodesic flow does not explode in {a = 0}, we are left with the latter possibility. It 
cannot lead to explosion either, since the geodesic flow needs a traveling time bounded 
away from to travel from W U W to V. > 

Note that the volume measure Vol£- of the compact manifold £ is an invariant finite measure 
for the modified diffusion. 



5.4. Crossing times and escape rate of 9-diffusions. Fix a reference point $ref £ 
and set P(-) = P($ref , ■)• L^^ ^^ emphasize that P is a two points function, so it is easy to 
pass from P($ref , ^) to P($o i^), ot the other way round, using the triangle inequality, for 
any $o G 
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Given an increasing sequence (-R„)n^i of positive reals, set Tq = and associate to each 

the exit time r„ from the box B^"^' := JD ^ Rn}- 

It takes the diffusion an amount of proper time (r„ — r„_i) to go from the box B^"'^^^ to the 
box B^"'\ The strategy in |HQ10| is to estimate P$(r„ — r„_i ^ t„) for a suitably chosen 
deterministic sequence {tn}n^Q of increments of time. Set for n > 1: 

n 



k=l 



If one can show that 



(5.4) 5Z^*(^« ~ ^"-1 ^ in) < oo 

for a convenient choice of the sequences {Rn)n^i and (T„)„^i, then the Borel-Cantelli lemma 
tells us that the diffusion does not exit 5^"^ before time T„ , for n large enough, preventing 
explosion. Following |HQ10| , we are going to consider the events 

SO as to be able to use our modified process run backwards from the fixed time T„ , when 
estimating the probability that the process crosses from 5^"^^) to 5^"^ not too fast. Lemma 
2.1 of |HQ10| (an application of the Borel-Cantelli lemma) justifies that considering these 
events leads to the same non-explosion conclusion as (15. 4p . We recall it here for the reader's 
convenience. 

Lemma 18 ( IHQIO Q. Fix $ e OM . // Yl P#(-E„) < oo , then there exists ¥^- almost- surely 

6 such that Tn ^ Tn — S, for all n ^ 1. 

We shall use the results of Sections I5.L1I and 15.31 to prove the fundamental estimate of 
Proposition [19] below. Given any compact subset B of OM, denote by P^ the law of the 
relativistic diffusion in OM started under the uniform probability in B: 



^^(■) = ^r47^ /p$(-)Vol(c?$). 



Similarly, and given any compact subset A oi £ , write Qa for the law of the modified 
0-diffusion in £ started under the uniform probability in A. 

PrDPOSITION 19. There exists a constant C such that we have for any n ^ 1: 



(r„. - 1 - 4t„) 

OZ KJn tn 



, . VOl(5W) TnJQn/tn 

where 0„ denotes the supremum of over the box \T> ^ /?,„ + l}. 

The proof mimics Takeda's original proof, as adapted by Hsu and Qin in |HQ10 , with the 
noticeable difference that we are working with a non-symmetric, non-elliptic diffusion. 
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Proof - We start by embedding the box B^"'^ into the set E^""^ constructed in Section [5731 
with A = Rn and e = | , say. From now on we work on the path space over S^"'^ and use 
the coordinate process X, whose filtration is denoted by (-Fs)s^o- We still denote by r^ 
the exit time from (the image in E^'^'^ of) B^^'> ; the event 

belongs to Jv^ . As explained above in Section 15.21 the proof has two main ingredients, 
the first of which is Inequality ( 15.51) below, where Q£-(n) denotes the distribution of the 
modified 0-diffusion in S^'^\ with generator Q given in (15. 3p . 

As the B-diffusion and the modified 0-diffusion have the same law before the stopping 
time r„, we have P^(„)(E„) = QB(n)(-E„) < 2Q£(„)(E„), and so 

Vdt ( R^")"! 

(5.5) Pb(i) (^n) ^ 2 ^^^^^ Q^(") {E^) , 

by the obvious inequality P5(i)(_E„) ^ yoUbW) JP'g(")(-^n)- The second ingredient involves 
the Lyons-Zheng decomposition of T>{Xs) under Q^cn). As T> is not a priori sufficiently 
regular to use Ito's formula, we apply it to its smooth approximation F constructed in 
Proposition [T2] (with R = Rn and t] = i). As the process (Xt'„_s)o<;s^t„ is under Q£-(n) 

d 

a homogeneous diffusion process with generator Q* = —Hq + ^"^Vji^aQ Vj) , it follows 

i=i _ 

from Ito's formula that there exists two martingales (Ms)o^s<gT„ and {Ms)o<:s^t„, with 
respect to the forward and backward filtrations of the process respectively, such that 

F{Xs) = F{Xo) + Ms+ f gF{Xr) dr, 

Jo 

F{X,) = F(Xr„_(^„_,)) = FiXrJ + Mt^s + j " Q*F{X,) dr, 
with 



(5.6) 



{M)s = V / aiXr) e(X,) \VjF\\Xr)dr ^ 4 0„s, 

„-i ■''0 



d 



^>s = E / «(^T„-r)e(Xr„_,) \V,F\\XT,^.r)dr ^ 4e„s. 



Setting M^ := Mt„-s and noting that Q — Q* = 2Hq, we thus have 

(5.7) d{F{Xs)) = d[ ^'^^' ) + HoF{X,) ds, 

with a controlled drift term \HqF\ ^ 2 , by Proposition [T21 By construction, we have 

sup |F(X._,+,)-F(X,_J| ^r„-l 

on the event E.^ , where X hits the set {F ^ Rn — \} in the time interval [t„_i, t„_i -|-t„]. 
To control the Q£-(n) -probability of En-, we use Hsu and Qin's trick. Cut the interval 

[0,T„] = U [(^ ~ l)^n,^^n] into £„ := T„/t„ sub-intervals of length t„ (to lighten the 

fc=i 
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notations, we shall neglect the fact that £„ may not be an integer ; this fact causes no 
trouble but notational), and write on each event {{k ~ l)t„ < r„_i < A;t„} 

This simple remark shows that the event < sup |F(Xt-„_;^+s) — F(Xt-„_J| ^ r„ — 1 > is 



0^ss£t„ 



included in one of the £„ events < sup \F{Xkt^^s) ~ F{^ktn)\ ^ ^'^2' f ' '^here 1 ^ /c ^ 

£n- By (15. 7p and the inequality \HqF\ ^ 2, the k^^ of these events is included in the union 
A^U Aj. , where 

Ak := I sup \Mkt,,+s - Mh„ I ^ ^ - 2t, 
Us£|s|s:i„ 

and 

I, := ( sup |MI,„+, - ii^^J ^ i:^ _ 2t„ 

Let W he a Brownian motion defined on some probability space (fi, J-", P). By (15. 6p we 
have 



h^.Mk) ^ 2P( sup IW^.I ^ --^'j ^ ^ V "1" ^^P f- ^""" ^- ^^"^' 
\o^s^t„ 4Ve„ y r„ - 1 - 4t„ V 32 6. t 



'n "n 



for some positive constant C; the same identity holds for Ak, using fl5.6p . Summing 
over A; and using Inequality (15. 5p yields the statement of the proposition since En C 

uiii(Afcuifc). > 

This key proposition being proved, it becomes easy to prove theorem [131 

Proof of Theorem [H - Taking i?„ = 2"+^ and t„ ^ 2"+^ in Proposition [TH so that 
T„ ^ 2"+^, we get for any n ^ 1 : 



Vdt f R'-")'! / O 
(5.8) P^d) (K) = Pij(i) (r„ - r„-i ^ t„, r„ ^ T„) ^ C ^^ )' \ -^ exp 

Vol ±3^ -' V tr, 



4" 



Specifying the choice of t„ by setting 

An-l 

tr, := min <^ 2"+\ 



1 + Iog+[0„VOL(5W)] 0„ 



the right hand side of (15. 8p is seen to be bounded above by a constant multiple of 2 "■, 
ensuring as a consequence the convergence of the series ^ P^(i)(ii^„). Indeed, we get 

from (15. 8p . with the above t„. 



,(E„) < C-VOL(i?(»)) J®- ^"g [Q-V"^(^^"^)] ,-41o4e.VOL(.(^))] < ^./2n_ 



(Ignoring the trivial case O = 0, we can suppose without loss of generality that we have 
G„ VoL(i?*^"'') ^ 3 for n large enough.) Note that the above choice of time increments t„ is 
simpler than Hsu and Qin's choice in |HQ10| ; there is in particular no need to introduce 
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their auxiliary function h{R) = log log i?, to get Grigor'yan's criterion, if the second 
upper bound of their section 3 is not used. 

n 

To conclude that the B-diffusion does not explode we need to check that T„ = Yl ^k 

fc=i 
increases to infinity. For the above choice of time increments tn, we have P^(i)-almost- 

surely, for n larger than some Hq, and for a positive universal constant c: 

n ( ,k_i 



T > \^ min < 2^=+^ - 

^~ £:. I ' ©2^+^+1 ( log^ [e2.+5+iVOL(i?2.+5)] + 1) 




(5.9) > c / min < 8 , . — -^ } dr 

^ ^ - ' ^ e, log [e, vol(5,)] ' 

Leaving aside the trivial case = and recalling that the map r i— )■ O,. = max 6 is 

non- decreasing, we can suppose without loss of generality that 0^ ^ 3. The divergence 
of the sequence (T„,) is then granted by the integral criterion 

dr = oo. 




Qr log [&r yOL{B, 

As Gj. increases, this condition is equivalent to 

n 

min < 8 , 1^ 

I «„ Ino- \f 

that is to 



> min < 8 , p ; -^ S = oo, 

^ 1 e„iog[e„voL(B„)] I 



E 



n 

oo. 



^ e,iog[e„voL(5„)] 

since the former holds obviously if an infinite number of terms were larger than 8. The 
previous condition is equivalent to Condition (15. 2p of Theorem [131 

Using Borel-Cantelli lemma under the form of Lemma [181 it follows that we have 



(5.10) P^d) sup V{^s) < '^ for any large enough n = 1 , 

so sup ^($s) < C)0, for all t > 0, since T„ increases to oo. Would a realization of the 

path $s explode by time t, its projection in M would provide a timelike path with an 
accumulation point (for it stays in the projection of a compact set by hypothesis (H)), 
contradicting the strong causality assumption on M. 

To prove that the same happens under any P$g , notice that since the non-explosion event 
E belongs to the invariant a-algebra, the function OM 9 $ i— )■ P<j,(-E) is ^e-harmonic, 
hence continuous, as Qq is hypoelliptic. It follows that since 






the probability P<j.(-E') must be equal to 1 for all $ G i?'-^^ But as the ball -B^^-* was 
arbitrarily chosen, P<j)(-E') is identically equal to 1 everywhere. > 
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5.5. Upper rate function. Using essentially the same reasoning as in Section 4 of |HQ10| , 
the above proof yields almost for free the upper rate function for the 0-diffusion given in 
corollary [I4l See also |Gri99] for related results. We keep the preceding notations. 

Proof of Corollary [TH - We follow the argument of |HQ10| , Section 4, making sure that 
it works here as well with our choice for t„, and without their auxiliary function log log. 
Suppose first G non-identically null and recall inequality (15. 9p . in which we can forget to 
take the minimum with 8, by Proposition [20] below. By fl5.10p . this yields the alniost- 
surely inequality 

sup D($,) ^ 2"+^ 

0^s^c/i-i{2")-<5 



that is 



sup V{(!?s)<52R, 

0^s^ch-i(K)-5 



for large enough R. Letting R = h{{t + 5)/c), this entails sup T'($s) ^ 32 /i((t + (5)/c), 

hence sup 'P($s) ^ 32 h{Ct), for large enough t. This shows the claim under the 

probability P^ci) , and then under P$q as well, by the same argument already used at the 
end of the proof of Theorem [T3i Finally, in the geodesic case (0 = O) , the same holds 
with T„ > c2" = c/i(2"). > 

5.6. Estimates of the volume of the sub-Riemannian boxes and application. Let 

us begin with a crude lower estimate of the volume of the boxes B^. based on the vertical 
expansion in the 5*00(1, (i)-fibre of OM, without taking into account the horizontal expansion 
which depends on the curvature of the base Lorentzian manifold M. We used this lower 
bound in the proof of Corollary [TU 

Proposition 20. We have liminf — —^ > d - 1 . 

Proof - Fix a relatively compact neighbourhood U of rrio in M, above which OU is trivi- 
alized in W X 5*00(1, d). Assume without loss of generality that $0 corresponds to {rriQ, 1). 
By the ball-box theorem (see e.g. [MQn02j ) , the box Br = {T> < r} contains a neigh- 
bourhood V X B{l,e) of $0 ) for some e > and for r larger than some fixed ri. Using 
this argument a finite number of times, together with the triangle inequality for V, we 
see that the box {V < r} contains any neighbourhood U x B{1, g) of $0, for any ^ > 0, 
provided r is large enough, say no less than ro = ro{U, g). Take g larger than the diameter 
of SO{d). 

We easily see that the boxes {V < r} dilate in the vertical directions Vi, . . . ,Vd with speed 
r, as r increases. So {V < r} contains the product of U by the ball of radius (r — ro) 
in 5*00(1,6?) for r large enough. This provides a lower bound on Vol({X' < r}) by some 
constant multiple of the volume of the hyperbolic ball of radius (r — ro), from which it 
follows that there exists some positive constant c such that log VoL(i?r) > {d — l) r + logc, 
for r large enough. > 

To close this work, we give a non-explosion criterion involving only the geometry of M, 
rather than the geometry of OM as it appears in Theorem [13] through the sub-Riemannian 
boxes Br- 
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Proposition 21. Fix $o = (mo,eo) e OM, and define the Sq>,^-radius p|(,(m) of any 
m e M as the infimum of the S^^^-length of C^ paths joining itlq to m. Define the S^^^-ball 
i?f^(r) of radius r as the set Bq,^(r) := |m G M | p^^^(m) ^ r}, and set 

V^^(r):=VoLM(5|^(r)). 

r/ien there exists a constant C such that we have for all r > 

logVoL(5,) <C + {d~l)r + logV^{Ce'). 

Note that the S'$,)-balls B^^{r) and their volume depend only on the choice of $o = 
(mo, eo) e OM and on the geometry of M. We noticed indeed in Section H] that the ^ey- 
length of a path in M started from m-o is the Euclidean length of its anti-development in 

(T^„M,eo). 

Proof - By the definitions in Sections H and IS.l.lj the b-distance of $o to any $ G OM 
is not larger than T'$g($), so B^. C B^{^Q;r), where B^ denotes the ball in OM of the 
b-metric. Vertically, that is to say in the frame tJ_^^($o) parallely transported along a 
minimizing curve 7, the maximal hyperbolic distance reached by the velocity component 
TJis of 7s is s, which is responsible for a maximal vertical volume (9(e^'^~^^'"). 
Having accelerated till reaching a maximal velocity 0[e^\ a minimizing curve in i?^($o; '") 
can perform a maximal horizontal displacement 0[e!'\ Hence we have the inclusions 

<(r)C7ro(i?^(<l>o;r))c<(0(e^)), 

and so Vol(5,.) < C e^'^-^> V^{Ce). > 

Applying Proposition [21] to the integral condition of Theorem [13] yields in the case of 

/"°° rdr 

a bounded the non-explosion criterion / ; ^r, — r = 00. Using the increasing 

y r + logy^(e^) 

character of the map (r h-)> V'^(e'')), discretizing and distinguishing whether or not there are 

infinitely many n such that log 1^*^(6") ^ n, we easily see that this condition is equivalent to 

the condition /°° ^^^ys^^^r-^ = 00. 

Corollary 22. Let {M.,g) be a strongly causal Lorentz manifold satisfying the complete- 

f°° rdr 

ness assumption (H) and the volume growth condition: / ■; ^. — - = 00. Then all 

J logV^^^e'') 

Q-diffusions with a bounded are stochastically complete. 

It is easy to see that this volume growth integral criterion does not depend on the choice 
of $0 £ OM. Contrary to Proposition [20l it relies on the horizontal expansion and not on 
the vertical expansion. This criterion does not apply to Godel universe, for which log ^'^(e''") 
is of order e*"; the non-explosion criterion of §3.21 covers the case of that spacetime. Corollary 
[22] applies for example to Lorentz manifolds which are topologically M}'^'^ and have a pseudo- 
metric g such that g,g^^, and the first order derivatives of g with respect to the canonical 
coordinates are bounded, since then logl^'^(e'') is of order r, as is the case in Minkowski 
spacetime. 
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